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ABSTRACT 

In  this  paper  we  develop  an  adaptive  beamsummer  for  di¬ 
rection  of  arrival  (DOA)  estimation  of  slow  fading  Raleigh 
signals,  using  a  new  design  approach,  termed  “Adaptive 
Detection/Estimation  for  sPecific  Tasks”  (ADEPT).  For 
DOA  estimation,  the  ADEPT  method  yields  a  weight 
adaptation  criterion  which  is  optimized  for  those  weights 
that  minimize  the  Cramer-Rao  (CR)  lower  bound  on 
achievable  mean-square-error  of  any  unbiased  DOA  esti¬ 
mator  constructed  on  the  beamsummer  outputs.  Simu¬ 
lation  results  are  provided  which  show  that  the  ADEPT 
DOA  beamsummer  yields  DOA  estimates  whose  mean- 
squared  error  (MSE)  approaches  that  of  the  more  complex 
maximum  likelihood  implementation. 


w' 


Figure  1:  Beamsummer  Model. 


1.  PROBLEM  STATEMENT 

Figure  1  shows  an  m-element  linear  beamsumming  ar¬ 
ray  on  which  impinge  p  signals  arriving  from  directions 
8\,...,8p.  At  the  fc-th  time  instant  a  snapshot  Y_k  — 
[Yi*, . . . ,  Y^]t  of  the  m-sensors  is  acquired,  k  =  1, . . . ,  n. 
We  assume  the  following  model  for  slow  Rayleigh  fading  [1] 
in  the  received  signals:  for  each  signal  the  amplitudes  are 
perfectly  coherent  over  the  array  for  any  particular  snap¬ 
shot  but  vary  randomly  for  different  snapshots.  In  addition 

we  assume  that  the  ambient  noise  N_k  —  [Nk . N„]T  is 

spatially  incoherent  and  temporally  uncorrelated.  Under 
the  above  assumptions  we  have  the  model  for  the  beam- 
summer  output  xk : 

xk  =  WHV(0)Sk  +  WHJV*  (1) 

where  5  =  [si , . . . ,  sp]T  is  a  vector  of  complex  signal  ampli¬ 
tudes,  and  V  is  an  m  x  p  matrix  of  steering  vectors.  With 
this  model  the  mean  and  covariance  of  the  k>h  snapshot 
are: 


px=E[xk]  =  WHV^ 

<rx  =  cov(x)  =  WwV£sVhW  +  WEjv  W  (2) 
where  E n  —  cr^I  is  the  covariance  of  the  array  noise. 

2.  OPTIMAL  BEAMSUMMING  FOR 

DOA-ESTIMATION 

Here  we  review  the  ADEPT-DOA  method  of  optimal 
weight  design  for  DOA  estimation  using  Fisher  informa¬ 
tion  and  the  Cramer-Rao  (CR)  bound  [2]. 


Let  <j>  be  an  unknown  parameter  of  interest.  For  fixed 
beamsummer  weights  W  let  the  N  snapshots  of  the  beam- 
summer  outputs  x  —  [a:1 , . . . ,  xw]r  have  joint  probability 
density  function  /*(x;  <j>)  parameterized  by  the  unknown 
scalar  parameter  <t>.  Define  the  Fisher  information: 


=  E$ 


^  In  4(x ;  <p) In  fx(x:  <j>) 


(3) 


Note  that  the  Fisher  information  is  a  function  of  both  <t> 
and  W.  Under  appropriate  regularity  conditions,  any  un¬ 
biased  estimator  <t>  =  <f>(x)  of  4>  has  variance  which  satisfies 
the  lower  bound  var[^]  >  -=S- .  The  methodology  of  optimal 
design  of  experiments  attempts  to  maximize  the  Fisher  in¬ 
formation  F$  =  FV(W)  over  W  in  order  to  provide  data  x 
which  allows  the  most  accurate  parameter  estimation  per¬ 
formance.  For  a  single  signal,  the  Fisher  information  for 
the  signal  direction  of  arrival  <j>  —  0  can  be  reduced  to  the 
form  [2]: 


where  u  =  D/X,  D  is  the  interelement  spacing,  A  is  the 
signal  wavelength,  and 


V  = 


■  21-m-l 

e3  2 


u  sin  0 


m 


(5) 
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=  [V£"^+lm^]2W, 
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2  m|/is|2 
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Note  that  S2  is  monotone  in  n2CJ,  which  is  the  maximum 
signal  variance  to  noise  ratio  of  the  output  of  the  beam- 
summer.  The  maximum  a2  is  achieved  by  the  adapted 
Frost  beamsummer  [3].  Also,  7 2  is  the  maximum  squared 
signal  mean  to  noise  ratio  of  the  output  of  the  beamsum¬ 
mer.  The  maximum  is  achieved  by  the  adapted  Ap- 
plebaum  beamsummer  [4].  The  ADEPT-DO  A  method  for 
beamsumming  seeks  to  choose  the  weight  vector  W  which 
maximizes  the  Fisher  information.  By  maximizing  -Fe(W) 
with  respect  to  W  we  achieve  an  optimal  compromise  be¬ 
tween  a2  and  7 o  by  minimizing  the  CR  bound  on  unbiased 
DOA  estimator  variance  at  the  beamsummer  output^ 
Define  the  normalized  vectors  V  =  V/|||Vj|  and  V_p  = 
ilp/||Kp||-  In  [2]  we  show  that  the  optimal  weight  vector 
that  maximizes  F$  is  given  up  to  a  complex  scale  factor 
by: 


W 


,Pt 


Vp  +  aV, 


6l  >  7o/2 

61  <  T 1/2 


(12) 


W  f  given  in  (12)  is  unimplementable  since  it  requires 
knowledge  of  8  and  the  SNR-dependent  constant  a.  We 
develop  a  direct  implementation  method  in  [2]  which  shows 
that  maximization  of  Fe(W)  is  equivalent  to  maximization 
of  the  function  Fnn a fW): 

WW  )  =  (14) 

w  w 

where 

cov^|Y]W 

COV2  [y]B? COV_2  [y]  +  COV-2  [y]BlCOV2  [y] 

4m-cov'^[y]BiE[y]E[yH]Bfcov-i[ri 
WHGiW 

~  H  ~  ’ 

w  w 

1, 


w  = 

Gi  = 

G2  = 
A'i(W)  = 
Ka(  W)  = 


and  Bi  and  B2  are  the  following  m  x  m  complex  matrices: 


Bi. 

B2 


<HaSl=l,...  ,m 


(  (2/-ra-l)/2\ 

\  \J m(m2  -  l)/6  J 


(15) 

(16) 


Note  that  Fnn  jfWl  can  be  used  for  weight  adaptation 
by  substituting  sample  means  and  covariances  for  E|Y]  and 
cov[y]  and  maximizing  over  W.  An  iterative  algorithm  for 
this  maximization  is  given  below: 

Eigenmaximization  Algorithm 


where 

_+  1  1 26%  ~  To 

1  +  cto  V  +  To 


(13) 


The  following  interpretations  follow  directly  from  these 
results: 

•  The  optimal  weight  vector  consists  of  the  sum  of  the 
two  unity  normalized  orthogonal  components  V  and 
V_p.  It  is  easily  seen  that  W  =  V_  gives  a  beamsummer 
which  maximizes  mean  array  gain  in  the  direction 
of  the  signal.  On  the  other  hand,  W  =  V_p  gives  a 
beamsummer  which  is  a  “signal  nuller”  -  it  minimizes 
mean  array  gain  in  the  direction  of  the  signal. 

•  For  So  <  7o/2,  up  to  an  arbitrary  scale  factor, 

WoP.  =  i-e-  the  optimal  weight  vector  places 

a  null  exactly  at  the  signal  angle  9.  This  suggests 
that  for  this  case,  the  best  ^-estimator  performance  is 
achieved  by  trying  to  null  out  the  signal,  e.g.  by  min¬ 
imizing  beamsummer  response.  On  the  other  hand, 
when  the  signal  energy  is  concentrated  in  the  sig¬ 
nal  variance  (a2  !§>  To)  the  Fisher  information  indi¬ 
cates  that  the  optimal  weights  must  pass  some  signal 
energy  -  adding  or  subtracting  o,l£  from  V_p  accom¬ 
plishes  this. 


{  Initialize:  W(0)  =  init.  weights,  v  —  tol,  i  =  1. 
while  ||W(«  +  1)  —  W(t)||  >  v 
{ 


W(i  + 1) 
i  =  *  +  1 


max 

e— vector 


GiAT(W(*))  +  G2A'2(W(t)) 


} 

w  =  cov-Jfy]  w(i  + 1) 


This  algorithm  requires  finding  the  square  root  factors 

i  _I 

cov 2  \Y]  and  cov  2  [y],  e.g.  by  Cholesky  decomposition, 
and  finding  the  eigenvector  associated  with  the  largest 
eigenvalue  of  the  matrix  G 1  A'i  (W1  +  G-j  /i'jfW).  A  steep¬ 
est  descent  algorithm  can  also  be  derived  which  reduces 
the  computational  load  but  only  guarantees  convergence 
to  a  local  maximum  of  Fdoa(YL)  [5].  A  block  diagram  of 
the  DOA-optimal  beamsummer  is  shown  in  Figure  2. 

3.  NULL  PICKING  FOR  DOA  ESTIMATION 


We  use  al  very  simple  gated  null  picking  method  to  estimate 
the  DOAs  from  the  adaptive  weights  given  by  the  above 
eigenmaximization  algorithm.  We  begin  by  computing  the 
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Figure  2:  DOA-optimal  Beamsummer  Block  Diagram. 


discrete  fourier  transform  (DFT)  of  the  adapted  DOA- 
optimal  weight  vector  W,  which  yields  an  antenna  gain 
pattern  for  the  adapted  weights  [2,5].  Nulls  in  this  beam- 
pattern  represent  possible  signal  DOAs.  To  determine 
which  nulls  represent  signal  DOAs,  we  begin  by  imple¬ 
menting  a  signal  detection-optimal  (D-optmal)  beamsum¬ 
ming  method  which  is  developed  in  [2,5],  called  ADEPT- 
D.  The  D-optimal  beamsummer  maximizes  the  function 
Fn( W)  over  the  weight  vector  W: 


Fd(  W)  = 


WH  [Gi  Ar  ( W)  +  G2  R2  ( W)]  W  „„ 

w*w  1  ; 


where 


Gr 

=  E\Y]E\Yh], 

Gi 

-  COv[y]-(T^, 

Ki  (W) 

=  1  +  2(1  +  a2)  +  y 

A2(W) 

2 

=  OL  , 

2 

Wh[Gi]W  1 

7 

WH W  <r2N  ’ 

2 

Wh[G2]W  1 

Ct 

WH  W  <t2n  ' 

The  D-optimal  beamsummer  is  implemented  using  the  pre¬ 
viously  described  eigenmaximization  algorithm.  The  re¬ 
sulting  D-optimal  weight  vector  produces  peaks  in  the  D- 
optimal  antenna  gain  pattern  near  the  DOAs  of  the  incom¬ 
ing  signals.  We  search  the  DOA-optimal  beampattern  for 
nulls  which  lie  within  1  beamwidth  “gates”  of  these  peaks. 
The  deepest  null  within  each  of  these  gates  is  used  as  a 
DOA  estimate.  This  process  is  illustrated  in  Figure  3. 


4.  SIMULATION  RESULTS 

An  i.i.d.  Gaussian  sequence  of  zero  mean  noise  variates 
was  generated  with  variance  =  1,  as  was  an  i.i.d.  Gaus¬ 
sian  sequence  of  signed  variates  with  prescribed  mean  n , 
and  variance  cr,.  From  the  signal  and  noise  sequences, 
1000  realizations  of  groups  of  100  snapshots  of  the  output 
of  a  six  element  uniform  linear  array  were  synthesized. 
For  Simulation  I,  we  set  y\  =  60  (eq.  9),  a20  —  6  (eq. 


D-optimal 


Figure  3:  DOA  Estimation  via  Gated  Null  Picking  for 
a  Signal  Arriving  at  8  =  .25  Radians. 


Simulation 

7  o 

at 

e 

I 

60 

6 

.22  radians 

II 

6 

60 

.22  radians 

Table  1:  Simulation  Signal  Parameters 


11).  For  simulation  II  we  set  y2  =6,  a2  =  60.  A  sum¬ 
mary  of  the  signal  parameters  used  for  the  simulations 
is  given  in  Table  1.  Figure  4  shows  a  typical  realiza¬ 
tion  of  adaptive  beamsummer  beampatterns  for  Simula¬ 
tion  I  and  varying  numbers  of  snapshots  for  the  Apple- 
baum  and  DOA-optimal  beamsummers.  Also  shown  are 
the  ensemble  average  beampatterns.  The  vertical  dashed 
line  represents  the  true  DOA.  Note  the  sharpness  of  the 
null  in  the  DOA-optimal  beampatterns  near  the  true  sig¬ 
nal  DOA.  Figures  5-6  show  the  MSE  of  the  DOA-estimates 
for  the  Applebaum  peak  picking  and  DOA-optimal  gated 
null  picking  estimators.  For  comparison  purposes,  a  clair¬ 
voyant  maximum-likelihood  (ML)  DOA  estimator  was  im¬ 
plemented,  operating  on  the  raw  array  data.  This  estima¬ 
tor  was  implemented  by  maximizing  the  array  likelihood 
function  f(Y\nsi  Ls,  &n,9)  over  6  assuming  that  the  true 
values  for  g_,,  Es,  and  <t2n  are  known  without  error.  The 
MSE  of  DOA  estimates  obtained  using  the  ML  estima¬ 
tor  is  included  in  Figures  5-6.  Finally,  we  plot  the  min¬ 
imum  beamsummer  CR-bound  l/mnrwfdW)  on  MSE 
for  any  unbiased  estimator.  Note  that  the  clairvoyant  ML 
MSE  falls  below  the  beamsummer  CR  bound,  as  expected. 
However,  since  the  ML  MSE  is  only  slightly  below  the  CR 
bound,  only  a  small  loss  in  estimator  performance  is  in¬ 
curred  by  restricting  ourselves  to  a  beamsummer  based 
estimator  implementation.  Also,  note  that  the  MSE  in 
the  DOA-optimal  estimates  approaches  that  of  the  clair¬ 
voyant  ML  estimator  and  the  beamsummer  CR  bound  as 
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Figure  4:  Typical  Beampatterns  for  a  Single  Realiza¬ 
tion. 
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Figure  5:  MSE  In  DOA  Estimates:  Simulation  I. 


the  number  of  snapshots  n  becomes  large.  The  improved 
estimation  performance  of  the  DOA-optimal  beamsummer 
over  the  Applebaum  beamsummer  can  be  attributed  to  the 
DOA-optimal  beamsummer’s  being  based  upon  minimiz¬ 
ing  the  CR  bound. 

CONCLUSIONS 

We  have  introduced  a  new  method  of  beamformer  design, 
(ADEPT),  which  adapts  the  beamformer  weights  to  pro¬ 
vide  the  minimum  CR  bound  for  DOA  estimation  based 
on  the  beamsummer  outputs.  Our  simulations  show  that 
when  implemented  with  a  very  simple  post-beamsumming 
DOA  estimation  technique,  the  DOA-optimal  beamsum¬ 
mer  approaches  the  CR  bound  on  MSE  for  any  unbiased 
DOA  estimation  technique  implemented  at  the  beamsum¬ 
mer  output.  In  [2,5],  we  introduce  other  optimal  beam- 
summers  designed  using  the  ADEPT  method  for  the  pur¬ 
poses  of  signal  detection  and  estimation  of  constant  mod¬ 
ulus  signal  parameters.  More  extensive  simulation  results 
can  be  found  in  [5]  for  multiple  signals  and  wideband 
beamsummer  arrays. 
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